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and solved. It is shown that the solution of the new evolution equation can provide a 
selfconsistent matching of "soft" high energy phenomenology with "hard" QCD physics. 
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1 Introduction. 



In this paper we discuss the QCD evolution for the gluon density in a nucleus. The gluon 
density is the most important physical observable that governs the physics at high energy 
(low Bjorken x) in deep inelastic processes Dealing with nucleus we have to take 
into account the shadowing correction (SC) due to rescattering of the gluon inside the 
nucleus, which is the main point of interest in this paper. We show that SC can be treated 
theoretically in the framework of perturbative QCD (pQCD) and can be calculated using 
the information on the behavior of the gluon structure function for the nucleon. 

Our calculations were performed in the double log approximation ( DLA) of pQCD, or 
in other words in the CLAP evolution equations for the region of low x ( high energy 
). In the DLA we consider the kinematic region where asln(l/a;) ^(Q^/Qq) ~ 1 while 
Q;sln(l/a;) -C 1 and as ln((5^/QQ) <^ 1 as well as as ^ 1, where is the virtuality 
of the photon and as is the QCD coupling constant. In terms of the anomalous dimension 
( see the next section for details) it means that we restrict ourselves by considering only 
the leading term in the anomalous dimension, namely 7 being as <^ 7 <^ 1. 

The advantages of the DLA are (see Fig.|l| for notations): 

1) we can neglect the change in the distance between quark (gluon) and antiquark 
(gluon) rt during the passage of the qq ( GG) pair through the nucleus. This simplifies 
the derivation of all formulae for the SC and leads to an eikonal picture of classical 
propagation of the qq (GG) pair with high energy which receives independent kicks due 
to rescattering in a nucleus. 

2) the cross section of the qq (GG) pair with transverse separation can be expressed 
through the gluon deep inelastic structure function for a nucleon. 

We will discuss both points in the next section in more detail. However, it is worthwhile 
mentioning that the DLA allows us to obtain the simplest and closest expression for the 
SC in DIS with a nucleus. 

The main goal of the paper is (i) to present a study of the SC in the Glauber approach 
using available information on the gluon distribution in the nucleon based mainly on new 
experimental data from HERA and on the solution of the GLAP evolution equations [0] ; 
and (ii) to find the generalization of the Glauber approach which will give the theoretical 
basis for the selfconsistent description of the gluon structure function for nuclei. 

Considering the deep inelastic scattering (DIS), we have to answer two principal 
questions: (i) why and how the cross section of DIS with nucleus which is equal to 

2 

a{'j*A) = Aa{'j*p) at x ~ 1 changes its A-dependence and becomes (7(7* A) oc Aa at 
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Figure 1: The structure of the parton cascade in the Glauber formula. A denotes the 
nucleus, N - the nucleon, G*{Q^) - the virtual gluon and ^^{r'^) is the nucleon cross 
section. 
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X — ^ ; and (ii) why a{'y*A), which is proportional to Aa at small values of x and Q^, 
approaches to a{'y*A) = Aa{'y*p) at large values of even in the region of small x. The 
last statement is obvious from the intuitive physical picture because at large values of 
(t(7*p) is small and the virtual photon probes the number of nucleons in a nucleus. 

The first question has been answered in the framework of the parton model ( see refs. 
[00) and the answer is based on the space-time picture of DIS in the rest frame of the 
nucleus. Indeed, the incident electron penetrates the nucleus and radiates the virtual 
photon whose lifetime t^* oc We can recover three different kinematic regions: 

'^'y* ~ mx ^ Rnn, where Rnn is the characteristic distances between the nucleons 
of the nucleus. This virtual photon can be absorbed only by one nucleon and the total 
cross section is a (7* A) = Acr^Yp)- 

2. Ra > ^-y* = > Rnn, where Ra is the nucleus radius. In this kinematic 
region the virtual photon can interact with the group of nucleons. However, a{'~f*A) is 
still proportional to A since the number of nucleons in a group is much less than A. 

3. r^* = — > Ra. Here, before reaching the front surface of the nucleus, the 
virtual photon "decomposes" in the developed parton cascade which then interacts with 
the nucleus. It can be shown that the absorption cross section of the virtual photon 

2 

will now be proportional to the surface area of the nucleus (7(7* A) oc As, because the 
wee partons of the parton cascade are absorbed at the surface and do not penetrate into 
the centre of the nucleus. 

However, the above simple picture cannot help us to answer the second question. 
Indeed, we can use it to explain the A dependence of the initial partonic distributions but 
in the GLAP evolution equations the A dependence is factorized out and do not affect the 

evolution. Therefore, we have to change the evolution equations to incorporate the 
physical phenomenon which we have formulated as the second question. To recover the 
physical origin of the new evolution equations in nucleus let us consider the oversimplified 
structure of the parton cascade: the virtual photon decays only in quark - antiquark pair. 
In this case the cross section can be written in the form: 

a(7*,A) oc 1^ dz I d^rt^{z,rt)(TA{zQo,r^)^*{z,rt) (1) 

where ca is the cross section for qq interaction with the nucleus, z is the fraction of energy 
of the photon (Qo) carried by quark and rt is the transverse separation between quark 
and antiquark. is the wave function of the virtual photon, which is known and is 
equal to Kf{art)[z'^ + {1 — z)"^], where Ki is the McDonald function and = Q'^z{l — z). 
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The main contribution in the eq. (0) comes from the region where art < 1. Expanding 
Ki and integrating over z(l — z) < < ^ we derive 

a 7*, A ^ 7^ L -^M^,r^) 2 
Q rf Q^rf 

For very small the cross section in QCD is small and proportional to r^. Such small 
cross section leads to cta oc A-^lnQ"^, since the probe with small cross section interacts 
with all nucleons in the nucleus. One can see that this is the first term of the GLAP 
evolution equations. With more complicated parton cascade we are able to reconstruct 
the GLAP evolution equations in a full. However, even at small r^, cr^ = Aaj\f only if the 
rescattering of qq - pair in a nucleus is small. The parameter which controls the value of 
the rescatterings is the number of collisions which is equal = 2crjv pRa, where p is the 
nucleon density in the nucleus. If z/ <^ 1 we can neglect the rescatterings, but if z/ ^ 1 
(Ta = 2ttR\ and does not depend on nucleon cross section. Therefore, we can trust the 
InQ^ contribution in eq. @) only for v = t'^pRa < 1 or for > pRa- Really, this 
condition depends on x too, since the nucleon cross section is a function of x, namely we 
will show that oc r?xG(x, -3-). 

Therefore, the lesson learned from this simple exercise is that we can trust the GLAP 
evolution for the nucleus structure function only for > Qq{x). It means that we have 
to solve the GLAP evolution equation with the initial condition on the line = Qq{x) 
or we have to change the evolution equations if we want to solve them in an usual way, 
namely, starting from the structure functions at = Qq- 

The previous attempts to attack this problem were related to the GLR equation |]T], ^ , 
in which the interaction ( recombination ) between two partons from different parton 
cascades was taken into account ( see refs.p, 0, |10[ [ll], |12] ). It was shown that the 



GLR equation is able to describe the main features of the experimental data on the deep 
inelastic scattering off nucleus [|14|, 0. The applications of the same ideas to description 



of other process such as J/\& production and Drell-Yan process |16] was also with 



reasonable success \l7,n 



However, the GLR equation was derived in the limited kinematic region where more 
complicated recombination processes have been proven to be negligible. Roughly speak- 
ing, we can trust the GLR equation only for small x and big values of Q^. 

In this paper we reanalyse the situation with the shadowing corrections in QCD for 
nucleus gluon structure function starting with the Glauber approach. The Glauber for- 
mula for the gluon structure function in the nucleus has been proven by A. Mueller in 
ref.|]T^ but this was remained unnoticed by the majority of experts because his paper 
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was devoted to quite different problem. In section 2 we rederive the Mueller formula and 
discuss the main properties of the Glauber approach in QCD. 

In section ^, we present the result of our calculations based on Glauber approach and 
point out which information on nuclear gluon distribution is needed in order to provide a 
reliable calculation. We adopt to following semiquantative way to study the mechanism 
and value of the SC: we used the GRV parameterization for the gluon structure function in 
the nucleon, which describes quite well the current experimental data. However, we would 
like rather to study the main properties of the SC than to provide reliable predictions for 
an experiment. 

In section ^ we consider the correction to the Glauber formula and discuss the general- 
ization of the Glauber approach for deep inelastic scattering off nucleus. Here, we discuss 
the GLR evolution equation which was the basis of all previous attempts to go beyond 
the Glauber approach and suggest and solve a more general evolution equation which is 
correct in the whole kinematic region without the limitations of the GLR equation. 

A summary and final discussion are presented is section ^. 

2 The Glauber approach in QCD . 

In this paper we will use the following notation (see Fig.|l| and Fig.^: 

denotes the virtuality of the gluon in deep inelastic scattering (DIS); 
m is the mass of the proton nucleus specified otherwise; 

Bjorken x = Xbj = ^ where ^/s = W is the cm. energy of the incoming gluon. 

kt denotes the transverse momentum of the quark or gluon; 

ft is the transverse separation between quark (gluon) and antiquark (gluon); 

bt is the impact parameter of the reaction which is the variable conjugated to the 
momentum transfer (g^); 

lit denotes the transverse momentum of the gluon attached to the quark - antiquark 
(gluon- gluon) pair. 

We will use the GLAP evolution equations for the parton densities in momentum 
space. For any function /(x), we define the moment f{uj) as 

fiuj)= fdxx-f{x). (3) 
Jo 
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Figure 2: Kinematics of diffractive production of a gluon pair. Q* is the virtual colorless 
probe of the gluon density. 
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Note that the moment variable uj is chosen such that the u; = moment measure 
the number of partons, and the lo = \ measure their momentum. An alternative moment 
variable = cu — 1 is often found in the literature. The x-distribution can be reconstructed 
by considering the inverse Mellin transform. For example, for the gluon distribution it 
reads 

xG{x, q2) = 1 / rf^ ^(^^ g2) ^ 

lUl JC 

where g{uj,Q'^) is the moment of the gluon distribution. The contour of integration C is 
taken to the right of all singularities. 

The GLAP evolution equations have the solutions of the form 

g{u:,Q')=g{u)e^^-'^^-Q\ (5) 

where 7(0;) denotes the anomalous dimension, which in the leading ln{l/xB) approxima- 
tion (LLA) of pQCD, is a function of a^/uj and can be presented as the following series 

■n oj Tc^ uj^ 

where C(3) is the Riemman zeta function, and is the number of colors. In DLA, only 
the first term in the above series is taken into account. 

The amplitude is normalized such that 

^ = ^l/(^,t)p, (7) 

atot = 47r/m/(s,0) , (8) 
and the scattering amplitude in ft^-space is defined as 

a{s,h) = ^ I d\e-''^^ '^ f{s,t = -q^). (9) 



27r 

In this representation 

atot = 2 y (fbt Im a{s,bt), (10) 
aei = I d%t \ a{sM)?- (11) 



The normalization of the nucleus wave function \E'A(ri, ...Ti^rA) is 

A 

J ^A(ri, ...ri, rA^Ain, ...n, ta) H d^^i = A , (12) 



i=l 
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where A is the number of nucleons in the nucleus, and the nucleon form factor in the 
nucleus is defined as 

A 

FA{qz,bt)= / dzie"'^''<iJA{zi,bt,r2,...ri,rA)'^*Aizubt,r2,...ri,rA)Y[d^ri. (13) 

i=2 

Throughout this paper we use the Gaussian parameterization for FA{qz,bt), namely 

FA{q.,h) = —^e ^ (14) 
where the mean radius of the nucleus R\ is equal to 

2 2 2 

Ra = '^R'WS ' (15) 

and RiYs is the size of the nucleus in the Wood-Saxon parameterization , which we 
choose R^Q = tqA^^^ with tq = 1.3/m. 



In nonrelativistic theory for the nucleus we can neglect the change of energy for the 
recoil nucleon. Its energy is Epi = m + g^/2m in the rest frame of the nucleus and 



2.1 Passage of the qq {GG) pair through the target 



The idea how to write the Glauber formula in QCD was originally formulated in Ref. |22 



and it has been carefully developed in Ref.|T9|. It is easier to explain the main idea 
considering the penetration of quark anti-quark pair, produced by the virtual photon, 
through the target. While the boson projectile is traversing the target, the distance rt 
between the quark and anti-quark can vary by amount Art oc Ra^ where E denotes the 
energy of the pair in the target rest frame and Ra is the size of the target (see Fig.|^). 

The quark transverse momentum is kt oc 1/rt. Therefore 

Art (X ^ < (16) 

and is valid if 

r^s > 2m R, (17) 
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Figure 3: The structure of parton cascade in the Glauber ( Mueller ) formula. 
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where s = 2mE. In terms of Bjorken x, the above condition looks as follows 



X « (18) 
2m R 

Therefore the transverse distance between quark and antiquark is a good degree of freedom 



[|T^[^p4[. As has been shown by A.Mueller, not only quark - antiquark pairs can be 
considered in such way. The propagation of a gluon through the target can be treated in 
a similar way as the interaction of gluon - gluon pair with definite transverse separation 
Tf with the target. It is easy to understand if we remember that virtual colorless graviton 
or Higgs boson is a probe of the gluon density. 

The total cross section of the absorption of gluon(G'*) with virtuality and Bjorken 
X can be written in the form: 

crfotiG*) = (19) 

f'dz [ ^ [^^f{Q\n,x,z)2{l - exp[-la{r^) S{b',)]} [^f{Q',n,x,z)r , 
Jo J 2tt J 2n 2 

where z is the fraction of energy which is carried by the gluon, \E'j^* is the wave function 
of the transverse polarized gluon and cr(r^) is the cross section of the interaction of the 
pair with transverse separation with the nucleon, and S{bt) is the profile function of a 
nucleus which we will specify later. 

The physical interpretation of Eq.(17) is very simple, if one notices that the factor 
in curly brackets is the total cross section for the GG pair with transverse separation 
passing through the nucleus 

rSC(rA = 9 / ^{1 - e" "^'^^^'-'^ } . (20) 

TT 



Indeed, the above formula is a solution of the s-channel unitarity relation: 

2Ima{s,bt) = \a{s,bt)\^ + G^sA) . (21) 

where a denotes the elastic amplitude for the GG pair with a transverse separation r^, 
and Gin is the contribution of all the inelastic processes. The inelastic cross section is 
equal to 

.... ^/<i'^*,G.„(..M^ /A (l-e-™<-^-) (22) 

We assume that the form of the final state is a uniform parton distribution that follows 
from the QCD evolution equations. Note that we neglect the contribution of all diffraction 
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dissociation processes to the inelastic final state (in particular to the "fan" diagrams which 
give an important contribution), as well as diffraction dissociation in the region of small 
masses, which cannot be presented as a decomposition of the GG wave function. We 
evaluate this input hypothesis in section H. 



In the language of Feynman diagrams, Eq.(pOD sums all diagrams of Fig.|^ in which 
the GG pair rescatters with the target, and exchanges "ladder" diagrams each of which 
corresponds to the gluon density. This sum has already been performed by Mueller [0, 
and we will only comment how one can obtain the result, without going into details. We 
borrow the presentation of the end of this subsection as well as the next two ones borrow 
Ref.|2|. 



The simplest way is to consider the inelastic cross section (see Refs.|jT9|, |T], ^), which 
has the direct interpretation through the parton wave function of the hadron (see Fig.|^), 
as all partons that are produced are on the mass shell. In leading Inl/x approximation 
we have two orderings in time: 

1. The time of emission of each "ladder" by the fast GG pair, which should obey the 
obvious ordering for n produced "ladders" (see Fig.|^) 

ti > t2 ... > U > U+i ... > tn] (23) 



2. Each additional "ladder" which should live for a shorter time than the previous one. 
This gives a second ordering (see Fig.^]): 

ti-t[ > t2-t'^ > ... > ti-t\ > ... > tn-t'^; (24) 

Each "ladder" in the leading log approximation is the same function of tj — t- which we 
denote as a{t — t', r^). This fact allows us to carry out the integration over ti and ti — t'^, 
which gives for n emitted cascades: 

= J^f / * ^(^1 - ^i) ^"(^1 - '^i' ) • (25) 

For the case of the GG pair, both last integrations are not logarithmic, and in the leading 
log approximation we can safely replace the above integral by 

c^n = A-2^''(t-t',r^), (26) 
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where t — t' is of the order of — due to uncertainty principle, where q\\ = —. One n\ 



<?ii 



is compensated by the number of possible diagrams, since the order of the "ladders" are 
not fixed. Therefore, the contribution of the nth "ladder" exchange gives 



1 



(27) 



Applying the AGK cutting rules[^ we reconstruct the total cross section which results 
in Eq.( 



2.2 a{ru qf) at qt = 0. 

The expression for a{rt, q^) at qt = 0. was first written down in Ref.[^ ( see Eq.(8) of this 
paper). It turns out that a can be expressed through the unintegrated parton function 
first introduced in the BFKL papers[^ and widely used in Ref. The relation of 



this function to the Feynman diagrams and the gluon density can be calculated using the 
following equation: 



l2 



as{Q^)xG{x,Q^) = / dl"^ as{lf) (p{x,l^) {2t 



Using the above equation we reproduce the results of Ref. [E^ , which reads 



c 

2^ 



Where = 

We evaluate this integral using Eqs.(^,([TO|) and ( p^T]) and integrate over the azimuthal 
angle. Introducing a new variable = rth, the integral can be written in the form: 



Evaluating the integral over ^ (see Ref. ||34|| 11.4.18) we have 



SCpas 2 [ du _^ , . , ^ .r^ 1 . ) r(7(^)) r(l - 7(0;)) 
= Nf^l Mc2^" ^^"^ {V{2-^{uW ^ ^ 

In the double log approximation of pQCD where 7(^0") <^ 1 the cross section for A^^c = 3 
reads: 

a{n,ql) = ^ vr^ ( xG^^^^{x, 1) ] (32) 
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This result coincides with the value of the cross section given in Refs. ^ (if we neglect 
the factor 4 in the argument of the gluon density). We checked that Eq.(2.16) of Ref. |]30| 
also leads to the same answer, unlike the value for a quoted in Ref.[16] (see Eq.(2.20)). 

In the case of the passage of the GG pair, Eq. (P^) reads (for Nc = 3) 

3as(4) / 4 \ 

a(r„ = 0) = -^n^r^ [xG%''''ix, A) j (33) 

Eqs (|32D and (^3|) are valid only on DLA of pQCD because we obtained these expressions 
considering lit ^ and neglecting the longitudinal part of the momentum li (Zj|| <C In). 
It is very essential to realize this limitation of our approach for all further applications. 

2.3 The bt dependence of the scattering amphtude with a nu- 
cleon 

To deal with the SC we need to know the amplitude not only at qt = 0, but at all values 
of momentum transfer, so that we can calculate the profile function of the amplitude in 
impact parameter space. The gluon density in a nucleon depends only weakly on qt in the 
GLAP evolution (see Ref. |l|] for details). Therefore, all the g^-dependence comes from the 
form factor of the G1G2 (qq) pair with the transverse separation and the form factor 
of the target nucleus. 

The form factor for G1G2 [qq) pair with the transverse separation rt is equal to 

where ki (fc-) denotes the momentum of the gluon 'i' before and after collision. Each of 
the wave functions is exponential, and a simple sum of different attachment of gluons 
lines gives 

Fg^gM) = e''^ - e^'^'-'') (35) 

We have absorbed the last factor in the expression for the cross section, while the 
transform gives qt dependence of Fg^g2- After integration over the azimuthal angle it has 
the form 

Fg.gM) = Jo{^-y) (36) 
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The target form factor cannot be treated theoretically in pQCD. In our problem it 
consists of the nuclear form factor and the nucleon distribution in the nucleus. For our 
purpose the phenomenological exponential parameterization for the nucleon form factor 
will suffice, namely 

F^{ql) = e-l^? (37) 



with the slope B = 10 Gev^^[^ , extracted from experimental data on hadron - hadron 
collisions, if we put the Pomeron slope a' = [^]. For nucleon distribution in a nucleus 
we use Eq.(0). Finally, the resulting bt distribution looks as follows 



where 



S{b',) =JS^ {{h - KY) Fa (g., K) ^ (38) 



1 , fbtrt\ _^?±i 



Considering R\ ^ B ^ one notice that it can be safely neglected the 6f-dependence 
in gluon nucleon interaction. Indeed, SNih) is a steep function of bt in comparison with 
FA{qz,bt)- It means that we can make the integral over 6J in Eq. (P^ putting b'^ = bt in 
F^Ailzjb't)- The result is 

SAibt) = F^iq,,bt) (40) 



2.4 The Qz dependence of nucleon density (gluon lifetime cut- 
off). 

To calculate the value of Qz we have to consider the process of diffractive dissociation 
pictured in Fig.(H). This process is the AGK cut with small multiplicity of produced 



hadron of the first diagram for the SC . We can find Qz from the obvious equation 

(g + qf = (41) 

which gives 



Qz = 7^ = —pr- + (42) 

2Qq 2Qq 



since Qz = Qo - ^ — ^ Qo 
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Now, we have to calculate through kt and the fraction of energy z (see Fig.|^). 



Using the technique of ref . ||3T[ , one obtains 



= ih + = ^ + = (43) 

Z (1 — -2) -2(1 — Z) 

To calculate z we have to consider the gluon structure function that enters in the value 
of awir^) (Eq.(|32|)). Indeed, using {k'^ -1)^ = (see Fig.(D) we have 

{k', - If = k'i - 2{k'^n +f = (44) 

since <^ k2, we have from Eq. (^) 

k'i = -k^ = 2 {k,n = 2 (^kok - {ko + ^)(/o + £)) 

P 

= -koT^ = -2zQox'm (45) 



which gives 



(46) 



2x'Qom 

Substituting z in the expression for we have = 2x'mQo. Finally 

g2 = (2;' + x)^m^ (47) 

where x' > x. 



Eg. (^71) has to be substituted in Eq. Ql^) to obtain the dependence of the nucleon 
distribution in a nucleus on q^. The physical meaning of this dependence is very 
simple. The life time of the virtual gluon is equal to 

^=Q-^ = 2^ 

If r is smaller then the size of the nucleus r < Ra the gluon cannot interact with all 
nucleons in a nucleus. The number of possible collisions is of the order of 

^ 2'mx ^ ^ 

where p is the density of nucleons in a nucleus. Eq.(|l^ and Eq . (^Tf) give us a practical 
way to introduce the finite life time of the virtual gluon in our calculation. 
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2.5 The wave function of virtual gluon. 



Here, we are going to discuss the last ingredient of the master formula of Eg. (pOD , namely, 
the wave function of initial gluon with virtuality Q^. Basically, it was done by A. Mueller 
in Ref.|T^, however, we will discuss it in this section to clarify the approximations that 
we have to do to get Mueller's answer. It is easier to discuss the colorless probe with 
virtuality that interacts with gluons than the virtual gluon. There is a number of 
such probes, for example, the graviton or Higgs meson. We do not need to specify what 
particle we use as a probe. We only need to write down the momentum structure for the 
meson-two gluon vertex. For a scalar particle such vertex has the following structure 



(50) 



where all notation is clear from Fig. 
meson to two gluons). 



see, for example Ref . ||33|| for the vertex of Higgs 



Q (Q,x) 

Tnmnry 




k, ((1-z) Q„,k,) 



Figure 4: The vertex of a virtual colorless probe ( Q* ) to two gluons. 

We intent to use a frame in which the nucleon with momentum p is essentially at rest 
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(p+ <ti Q+) and where 



Q = iQ+,Q-,Qt) = iQ+,-^,o). (51) 



To find the wave function we can use the technique developed in Ref.|3^. however, 
before doing so we need to specify the polarization of gluons that work in our case. The 
problem is that we know that the gluon which interacts with the target in Fig.(^) (^2 in 
our case) has the longitudinal polarization at high energy (small x ), while the second 
gluon indeed has been produced on mass shell and has only transversal polarization (see, 
for example, Ref. |T[] for a detailed discussion). It means that we cannot treat the first 
gluon as a real particle on the mass shell, where gluons have only transverse polarization. 
Therefore, before introducing the wave function of the probe we have to make use of the 
polarization vectors and the Weisszaker- Williams transform from longitudinal polarized 
gluon to transverse one. Indeed k'^iyC^^^ = 0, where e^^^ is the polarization vector of the 
H — W gluon. Therefore 

^1' ) ■ = ■ (52) 
Since only components along (5+ of vectors fcj'^j and k^j^^^ are big, we obtain 



or 



Q+e« = -^!L^ (53) 



Finally using Eq.(^) we can rewrite the vertex T^j_i, in the form 

z{l — z) 

where we sum over two possibilities: gluon 2 interacts with the target and gluon 1 interacts 
with the target. 

We anticipate that sufficiently small transverse distances(rf) will contribute to the 
processes {jsoft rt ^) . In this case all interactions between gluon 1 and 2 before the 
interaction with the target are small since they are of the order of as(^t ) and as{rt) <^ 1. 
Therefore we have to calculate only the contribution of two gluon state to the wave 



function and using the technique of Ref.[^] and Eq.(^J) for F^j, we have 



^,,u{z, kt) 
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^Jz{l-z) (Q2^(l-z) + fc|) 
Going to Tt- representation one obtains 

1 



z) = -V^V,Ko{ar) ^ (56) 

z(l - z) 



where = Q'^z{l — z) and Ko{ar) is the McDonald function. Making use of the properties 
of the McDonald functions we derive 



= , ^ < a K2{ar)^^ ^ \ (57) 

'z(l - z) { 



where fi, u = 1, 2. 



2.6 The Mueller formula. 

Now we have all ingredients that we need to derive the Glauber formula for the perco- 
lation of GG pair through the nucleus. Using the well known relationship between cross 
section and the deep inelastic structure function we can derive the following formula for 
xGa{x, Q^)JotN^ = 3{ see Ref.^) 

xGaM') = ^ f'dz [ ^|vI/(r„z)p2|l-e-^'^^"(-''-?)^(^?)| (58) 

TT^ Jo J n J n ^ ' 

where x' = ^.-ilq^^ ■ To specify the region of integration we have to put in Eg. ([58|) the wave 
function of Eq.(|57|) 

= 5]^(rt,z)^*(n,2;) = 

_L_ . j [a'K,_(ar) - '^^^ ' + (59) 

where represents the sum over polarization of the gluon. The main contribution in 
Eq.(p5|) comes from the region of small z {art <^ 1), where 

|v^(ri,z)p = A (60) 
zrt 
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which can be easily derived from the expansion of McDonald function at art -C 1. In all 
terms of the expansion, except the first one, the integral over z is convergent. 

The condition of a <^ 1 means that 



Introducing the new variable 



x' = ^ ^ , (62) 



2zQorfm 

instead of z, one sees that Eq.(|6T|) can be rewritten in the form 



x' > = XBj = X (63) 

2mQo 



Substituting Eqs.(|60D and (0) in Eg. (|58D , we derive the Mueller's formula 

The lower limit in integration comes from eq. (|6lD . 

It is easy to see that the first term in the expansion of Eq. (|6^ ) with respect to cr gives 
the GLAP equation in the region of small x. Using Gaussian parameterization for S{bt) 
( see Eq.(|l^)) we can take the integral over bt and obtain the answer (A^^^ = Nf = 3) 

xGa{x, Q') = ^ ^ /5 ^ {C + ln{Ka{x', r^)) + E,{Ka{x', r^))} (65) 



where C is the Euler constant and Ei is the exponential integral (see Ref.||3^ Eq. 5.7.11 
and 

[X , 



^g{x', rl) = -^x'G^^^^(x', -) (66) 



To understand the physical meaning of this equation it is instructive to write down 
the evolution equation for the gluon density. Indeed, 

d'xG{x,Q') 

(67) 



ain(l/a;)(91nQ2 



Ncas ^rOLAPf^ ^ 2 (-1)^ 1 ( 7,N,AasxG^^^^{x,Q') \'^' 
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The first term corresponds to the usual GLAP equations, while the second one takes into 
account the SC. It should be stressed that the term with A; = 1 (if treated as an equation 
[Q), is the same term that appears in the nonlinear GLR equation which sums the "fan" 
diagrams. This term has been calculated using quite a different technique |3^. The 
coefficient in front of the other terms reflects the fact that all correlations between the 
gluons have been neglected, despite the fact that gluons are uniformly distributed in the 
disc of radius Ra- 

Mueller's formula is not a nonlinear equation, it is the analogue of the Glauber formula 
for the interaction with a nucleus, which gives us the possibility to calculate the shadowing 
corrections using the solution of the GLAP evolution equation. Hence, this formula should 
be used as an input to obtain the complete effect of the SC, for the more complicated 
evolution equation, such as GLR[Q, or the generalized evolution equations (see Ref . |]40| ) . 

Calculating the subsequent iterations of the Mueller's formula gives us a practical way 
to estimate the value of the SC from more complicated Feynman diagrams, that have to 
be taken into account (diagrams such as in Fig.^). Certainly, the iteration of Mueller's 
formula is not the most efficient way to calculate the SC correction in the region of 
extremely small x, but it could give sufficiently reliable results for the HERA kinematic 
region. 

At first sight eq. ( p7D looks like the Operator Product Expansion. However, it should 
be stressed that the Mueller formula itself takes into account not only the high twist 
contributions to the deep inelastic structure function but also changes crucially the leading 
twist term since it contains the full integral over r^. 

2.7 An instructive example. 

Before discussing the numerical results of integration of the master equation (^) we will 
consider one example which gives the physical picture there originated. 

The first observation is that Mueller formula itself gives a natural cutoff for the large 
distance contribution and can be used as a model for the behavior of the gluon distribution 
in a target even at very small values of Q^. 

Doing all calculations in double log approximation of pQCD ( or in other words, using 
the GLAP evolution equation in the region of small x ), we have to assume that the 
anomalous dimension 7(0;) <^ 1. It means that we can confidently take the gluon density 
independent of the integral at the low limit = and only integrate over the factor 
in the expression of a^'-^ . Straightforward integration leads to the answer 

xGa{x,Q') = -j-^ / —RlQ' {C + ln{KG) + (1 + kg)E,{kg) + 1 - e-'^^} (68) 
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Figure 5: The interaction with nucleons that is not taken into account in the Glauber ( 
Mueller) formula. 
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where 

-c(x) = ^x'G«^^^(x',g^) (69) 

The last integration over x' has to be done numerically. 

One can see from Eq.(^) that the answer mostly depends on hq. For hq Eq.(^) 
gives: 

xGa{x. Q') = ^ Q')ln ( (70) 



TT 

12/ 



where Ql{x') is the solution of the equation = I- 

This result is very close to GLAP equations. The difference is due to the fact that one 
cannot substitute the gluon structure function at low limit of integration over for the 
first term of expansion for Eq. (p^D . This particular contribution is of the order of 1/'j{uj). 
It is necessary to improve Eq. (^) by adding a term 

TT Jx X' Jmax{Ql,Ql{x')} Q''^ 



71 Ix x' i ) ^uiax{Ql,QQ(x')}^ ^ ^ 

where Qq is the initial virtuality from which we start the GLAP evolution. Eq . ([fO|) 
illustrates the important point that the SC provides us a new scale in the evolution which 
crucially depends in x'. It means that the first correction for the evolution equations can 
be found, just introducing this scale, namely 

In Fig.^ are plotted the contours of k, for a nucleon target ( A = 1 in eq. ( |69D ) that give 
an idea in which kinematic region we expect big SC. 

The large SC corresponds to ^ 1 where Eq.(|68D gives 

xGa{x,Q') = -R\Q' ^{C + lniKc)) (73) 

TT Jx X 

Sum of Eqs. ([7^) and (|73D gives the answer 

xGAx, Q') = e{Q^- Qlix)) ^ ^ %x'G%^^^{x\ Q^) 

^ ^ n Jx x' JQl(,x') Q'^ 

+ e {Ql{x) - Q^) -R\Q' —{C + ln{Kc)) (74) 
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Figure 6: The contours of k for Nucleon, Ca and Au. 
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Eq.([74D is very approximate but Fig.(|^) shows that this equations describes the full 
formula quite well except the region where ~ Q\[x). The nice feature of this equation 
is that it illustrates in explicit way how the SC work. First, they provide the new scale 
for the transverse momentum inside the parton cascade [Qq[x)). It means that we expect 
the GLAP evolution only for distances < \IQq{x). Second, the SC generate the 
surface term ( the second in Eq.([74D ) which gives a oc J^l"^ as for normal hadron nucleus 
interaction. Of course, all these properties have been anticipated (see Refs.0, ||, ^) and 
Eq.([7^) gives the simplest example how they manifest themselves in the case of the deep 
inelastic scattering. 

Now, let us try to understand x' dependence of our integrand. In our formula it is 
implicitly assumed that each gluon can interact with all nucleons in a nucleus. It is not 
the case because each gluon lives a certain time which is tq = X/mx' where coefficient 
A is not quite well known, but A > 1/2. It means that the gluon can interact only with 
T = R\ptg nucleons unlike T becomes smaller than A. Substituting T instead of A we 
have the following answer for Eq.(^) 

3as f^dx'^^^,^ fQ^ dQ'^ ,r^GLAPrf n2\ 



7T Jx X' JQoi^') Q 

at x' > xa = X/mRA we have no log integration over x'. Therefore, basically, the GLAP 
equation can be reduced to the form 

3asA r^A dx' fQ' dQ'^ ,r.GLAP, , 



TT 



Q') = ^ — I ,^ Q^) (76) 

Ix X JQg(x') 



One can see in Eq.(|70D that a new cutoff in x appears in the GLAP equation which has 



been discussed two decades ago in Ref. |42, 43 . 



2.8 Theory status of the Mueller formula. 

In this section we shall recall the main assumptions that have been made to obtain the 
Mueller formula. 

1. The gluon energy (x) should be high (small) enough to satisfy Eqs.(|l^) and 
asln{l/x) < 1. The last condition means that we have to assume the leading ln{l/x) 
approximation of perturbative QCD for the nucleon gluon structure function. 
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2. The GLAP evolution equations hold in the region of small x or, in other words, 
as/n(l/rj) < 1. One of the lessons from HERA data is the fact that the GLAP evolution 
can describe the experimental data. 

3. Only the fastest partons {GG pairs) interact with the target and there are no 
correlations (interaction) between partons from the different parton cascades (see Fig.|^). 

4. There are no correlations between different nucleons in a nucleus. 

5. The average for GG pair-nucleon interaction is much smaller than Ra- 

The first assumption allows us to treat successive rescatterings as independent and 
simplifies all formulas reducing the problem to an eikonal picture of the classical propaga- 
tion of a relativistic particle with high energy {E ^ where /z is the scattering radius 
in the nuclear matter) through the nucleus. The second one simplifies calculations but 



we can consider the BFKL evolution instead of GLAP one. The third assumption is 
artifact of the eikonal approach and we shall discuss it in the following sections. The last 
two are usual assumptions to treat nucleus scattering. We have used the specific Gaussian 
parameterization for bt dependence. Also, one can easily generalize our formula in more 
general C3iSG, clS Wood-Saxon parameterization ||2T[|. 



3 pQCD calculations from the Mueller formula. 



We use the GRV parameterization for the nucleon gluon distribution, which de- 
scribes all available experimental data quite well, including recent HERA data at low 
xBj- Moreover, GRV is suited for our purpose because (i) the initial virtuality for the 
GLAP evolution is small {Qq ~ 0.256*6^^) and we can discuss the contribution of the 
large distances in MP having some support from experimental data; (ii) in this parame- 
terization the most essential contribution comes from the region where OglnQ'^ ^ 1 and 
aslnl/xBj ~ 1- This allows the use of the double leading log approximation of pQCD, 



where the MP is proven [^]. It should be also stressed here, that we look at the GRV 
parameterization as a solution of the GLAP evolution equations, disregarding how much 
of the SC has been taken into account in this parameterization in the form of the initial 
gluon distribution. 

However, in spite of the fact that the GLAP evolution in the GRV parameterization 
starts from very low virtuality ( Qq ~ 0.25GeV^) it turns out that the DLA still does 
not work quite well in the accessible kinematic region (Q^ > lGeV'^,x > 10~^). On 
the other hand, our master equation (see Eq.^lD is proven in DLA. Willing to develop a 
realistic approach in the region of not ultra small x {x > 10~^) we have to change our 
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master equation ([6^). We suggest to integrate Eg. ([67|) and substitute the small x kernel 
for the full GLAP kernel in the first term of the r.h.s. This procedure gives 

xGa{x,Q^) = xGA{x,Q^){eq. (64)) + AxG%'^^{x,Q^) 

The above equation includes also AxG%^^{x,Qq) as the initial condition for the gluon 
distribution and gives AxG%^'^ {x,Q^) as the first term of the expansion with respect to 
Kg- Therefore, this equation is an attempt to include the full expression for the anomalous 
dimension for the scattering off each nucleon, while we use the DLA to take into account 
all SC. Our hope, which we will confirm by numerical calculation, is that the SC are small 
enough for x > 10~^ and we can be not so careful in the accuracy of their calculation in 
this kinematic region. Going to smaller x, the DLA becomes better and eq. (^) tends to 



our master equation (64). 



3.1 The gluon structure function for nucleon. 

In this subsection we are going to check how eq. ( |77D describes the gluon structure function 
for a nucleon, which is our main ingredient in the Mueller formula. We calculate first the 
ratio 

J, ^ xG^{x,Q'){eq. (77)) 

xG%^^{x,Q^) ' ^""^ 

for A = 1, which is shown in Fig.^ From this ratio we can see the general behavior of 
the SC as a function of ln{l/x) and Q^. When compared to the GRV gluon distribution, 
the xG^ distribution presents a suppression which increases with ln{l/x) and decreases 
with the virtuality Q^. In the region of the HERA data, 3 < ln{l/x) < 10, and > 
2 Gel^^0, the SC are not bigger than 15 %. The SC give a contribution bigger than 20 % 
only at very small value of x, where we have no experimental data. 

In the semiclassical approach (see |I[]), the nucleon structure function is supposed to 
have and x dependence as 

xG^{x,Q') oc {Q2|<7> |i|<-> . (79) 

X 

We can calculate both exponents using the definitions 

dHxG^{x,Q')) 
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Figure 8: The SC for nucleon (A—1) as a function of ln{l/x) and , where ratio Ri 
compares xG^ with xG (GRV) distribution. 
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< 7 > 



d\n{xG^{x,Q^)) 



^1) 



The eq.(Bn) gives the average value of the effective power < u > of the gluon distri- 



bution, xG{x, 



oc X 



-<UJ> 



which is suitable to study the small x behavior of the gluon 



distributions. Fig.|^ shows the calculation of < u > the nucleon distribution calculated 
using eq.(^) and for GRV gluon distribution, both as functions of ln{l/x) for different 
values of Q^. From the figure, we can see that the effective powers of xG^{A = 1) and 
xG{GRV) have the same general behavior in the small x limit but the nucleon distribution 
is slightly suppressed. We calculate also, in the same kinematical region, the exponent 
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Figure 9: The effective power < u > calculated for xG^ {A = 1) and the GRV distribution. 



< 7 >, given by eq (PT]). This is the average value of the anomalous dimension, which 
describes the effective dependence of the distribution in variable. Figs.|l^ shows < 7 > 
for the nucleon and GRV distributions, indicating that the dependence is slightly 
soften by the SC. 

Comparing figures ^ and IC , we can conclude that even these more detailed charac- 
teristic of the gluon structure function have not been seriously affected by the SC in the 
nucleon case. 

We also use the CLAP evolution equations to predict the value of the deep inelastic 
structure function F2 from the xG^ gluon distribution. Summing the CLAP evolution 
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Figure 10: The effective power of Q"^ dependence calculated for xG^{ A = 1) and the GRV 
distribution. 



equations for each quark flavor, the function F2 may be written ^ 

- ^Ee;/; f i'- . (1 - (8.) 

where the sea quark distributions have been neglected in comparison with the gluon 
distribution. Fig.^ shows the prediction for F2 from xG^ and from the GRV distribution, 
compared with experimental data. As we can see, the magnitude of the suppression due 
to the SC is less than 10% in the region of the HERA data and this suppression is smaller 
than the experimental error. 



From the above results we can conclude that eq. ([77|) gives a good description for the 
gluon structure function for nucleon and describes the available experimental data. The 
MF provides a good description of the SC and can be taken as a correct first approximation 
in the approach to the nucleus case. 



3.2 The gluon structure function for nucleus. 

In the framework of perturbative approach it is only possible to calculate the behavior of 
the gluon distribution at small distances. The initial gluon distribution should be taken 
from the experiment. Actually the initial virtuality Ql should be big enough to guarantee 
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that we are dealing with the leading twist contribution. Our main assumption is that 
we start the QCD evolution with a small value of Qq considering that the MF is a good 
model for high twist contributions in DIS off nucleus. 

The scale of the SC governs by the value of ka, namely they are big for > 1 and 
small for /t^ < 1. Fig. 6 shows the plot of = 1 for different nuclei. One can see that the 
SC should be essential for heavy nuclei starting from Ca at the accessible experimentally 
kinematic region. 

Now we extend the definition of Ri for the nucleus case 

' AxG'G«^(x,Q2) ' ^^^^ 

where the numerator is calculated using eg. ([77|) . Figure |T^ shows the results for the 
calculations of -Ri as a function of the variables ln{l/x), InQ'^ and A^/^. Fig.[T^a presents 
the ratio Ri for two different values of and for different nuclei. The suppression due to 
the SC increases with ln{l/x) and is much bigger than for the nucleon case. For A = 40 
(Ca) and = lOGeV'^, the suppression varies from 4% for ln{l/x) = 3 to 25 % for 
ln{l/x) = 10. For A = 197 (Au) the suppression is still bigger, going from 6% to 35% in 
the same kinematic region. Fig. |lpD shows the same ratio for different values of for the 
gold. The suppression decreases with Q^. Figs. and d show the Ri ratio as a function 
of A^^^ and x for a fixed value of Q^. As expected, the SC increases with A. An interesting 
feature of this figure is the fact that the curves tend to straight lines as x increases. It 
occurs because, as x grows, the structure function xG{GRV) becomes smaller, and the 
correction term of ( [77D proportional to n dominates. Since n is proportional to A^^^, the 
curves behave as straight lines. The decrease of suppression with is illustrated in more 
detail in Figs. |l2|e and f which presents -Ri as a function of In for different values of x 
for Ca and Au, respectively. The effect is pronounced for small and x and diminishes 
as InQ"^ increases. 

This picture ( Fig.|l^ ) shows also that the gluon structure function is far away from 
the asymptotic one. The asymptotic behavior i?i — >■ 1 ( see Figs.|T2|e and f ) occurs only 
at very high value of as well as in the GLR approach ( see ref. ). The asymptotic 
A-dependence ( -Ri oc A's) ) has not been seen in the accessible kinematic range of 
and X ( see Figs. |^c and d and Table 1 ). This result also has been predicted in the GLR 
approach [^]. We want also to mention that parameterization Ri = Rin A~°' does not fit 
the result of calculations quite well for IGev"^ < < 20GeV'^ and 10"^ < a; < 10~^. 
For X ~ 10~^ the parameterization Ri = Rin — R' As with parameters Rin and R' for 
each value of Q^, works much better reflecting that only the first correction to the Born 
term is essential in the Mueller formula. 
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Figure 12: Ri as a function ofln{l/x), InQ^ and A^/^: a) Ri as a function ofln{l/x) 
for different nucleus and different values of A; b) Ri as a function of InQ"^ for different 
values of xb for Au; c) and d) Ri as a function of A^l^ for different Q^; e) and f) Ri 
dependence on for Ca and Au. 
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Table 1: Values of Rin and a for parameterization Ri = Rin A 







Q' = 


IGeV^ 


Q' = 


lOGeV^ 


X 




RlN 


a 


RlN 


a 


10- 


-2 


0.94 


0.0416 


0.98 


0.014 


10- 


-3 


0.92 


0.0616 


0.94 


0.034 


10- 


-4 


0.88 


0.094 


0.92 


0.0563 


10- 


-5 


0.8 


0.145 


0.86 


0.093 



We extend also the calculation of the exponents < a; > and < 7 > of the semiclassical 
approach for the nuclear case. We calculate the effective power of the nuclear gluon 
distribution < a; > using the expression 



< a; > 



d\n{xG'^{x,Q^j) 
d\n{l/x) 



FigJT^ shows the results as functions of ln{l/x) for different values of and different 
nucleus. The SC decreases the effective power of the nuclear distribution, giving rise to a 
flattening of the distribution in the small x region. 
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Figure 13: < u; > for different values of and A. 
It is also interesting to notice that at small values of Q^, the effective power tends to 
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be rather small, even in the nucleon case, at very small x. However it should be stressed 
that the effective power remains bigger than the intercept of the so called "soft" Pomeron 
even in the case of a sufficiently heavy nucleus (Au), for Q'^ > IGeV"^. Nowadays, 
many parameterizations with matching of "soft" and "hard" Pomeron have appeared 
triggered by new HERA data on diffraction dissociation These parameterization 



used Pomeron-like behavior namely, xG{x,Q^) oc a;~^('3 ). However, if the Pomeron is a 
Regge pole, u cannot depend on Q^, and the only reasonable explanation is to describe 
a;(Q^) as the result of the SC. Looking at Fig.|l3|we can claim the SC from the MF cannot 
provide sufficiently strong SC to reduce the value of u to 0.08, a typical value for the soft 
Pomeron [^, at least for > IGeV'^. 



The calculation of the effective value of the anomalous dimension 7 may help us to 
estimate what distances work in the SC corrections. This effective exponent is given by 



< 7 > 



d\n{xG^{x,Q^)) 



(85) 
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Figure 14: < 7 > for different and A. 



Fig.|T^ shows the results as functions of ln{l/x) for different values of and for two 
nuclei. We see that the values of 7 at ln{l/x) < 5, for both Ca and Au, is very close to 
the results for GRV and for nucleon case. At smaller values of x, the anomalous dimension 
presents a sizeable reduction, which increases with A. For ln{l/x) > 15, < 7 > tends to 
zero unlike in the CLAP evolution equations ( see Fig. 10 for the GRV parameterization). 
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Analysing the dependence, we see that < 7 > is bigger than 1 only for Q'^ = 0.5 GeV^. 
For = 1.0 GeV^, the anomalous dimension is close to 1/2, and for > 5.0 GeV^ it is 
always smaller than 1/2. 

Using semiclassical approach, we see that 



(86) 

and if 7 > 1, the integral over rt in the master equation (^) becomes divergent, concen- 
trating at small distances. 

If 1 > 7 > 1/2, only the first SC term, namely, the second term in expansion of 
the master equation, is concentrated at small distances, while higher order SC are still 
sensitive to small rt behavior. Fig.O shows that this situation occurs for > IGeV'^, 



and even for = IGeV^ at very small values of x. We will return to discussion of these 
properties of the anomalous dimension behavior in the next section. 

In subsection 2.4 we have discussed that the virtual gluon can interact with the target 
only during the finite time r (see eq. (^) ) undergoing pr < p Ra collisions. In the 



framework of the Glauber approach the easiest way to take into account the finite life 
time of the gluon is to include in our calculation the longitudinal part of the transferred 
momentum ( ) to a nucleon during the collision ( see, for example, Ref. |^8|] ). Using 
eq. (|I0| ) and eq. (|T7|), we can obtain: 

xGa{x,Q') = AxGm{x,Q') - (' I ^^L{q^)x'GM{x',Q") (87) 

TT Jx J X LJ 

71 J X X J Tt 

where ^ 

L{q,) = ^-^mHx + x'? _ 



Fig. 15 shows the result of our calculations. Comparing Fig. 12 with this picture. 



one can see that the finite life time of the virtual gluon affects the behavior of the gluon 
structure function only at sufficiently large x {x > 10~^ ) diminishing the value of the SC 
in this kinematic region. More interested in the small x behavior of the gluon structure 
function in nuclei we neglect the finite life time of a gluon through all calculations below. 
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Figure 15: Ri for CA and Au with gluon life time cutoff. 
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3.3 The behavior of o^g. 



We calculated the ratio 





ainiaino2(2;G'A(a;, Q )) 









which is equal to 1 in the case of the GLAP equation. In some sense it characterizes 
the value of the effective QCD coupling constant in a nucleus. This constant becomes 
smaller due to SC and in a restricted kinematic region we can treat the parton cascade 
as a parton system evolving according to the linear GLAP equations but with smaller 



coupling constant. The result is shown in Fig.|T6|. The physical meaning of the weakening 
of the coupling constant is very simple. At low-x, the gluon density increases and in a 
typical parton cascade cell there are many partons with different colors. The color average 
gives smaller mean color charge at small gluon densities. One manifestation of the parton 
interaction weakening is the bigger value of A QCD for nuclei at small x. From Fig.|16 
one can see that the GLAP evolution equation can be used to describe the evolution of 
the gluon structure function for the nucleon in HERA kinematic region ( x > 10~^ ) but 
only for sufficiently large photon virtualities { > 5 GeV^ ) . At smaller values of 
we have to develop more sophisticated approach for the adequate treatment. In the case 
of nuclei, the GLAP evolution equations with screened ( reduced ) value of as could be 
considered only as the first rough approximation even for large {Q^ > 5GeV^ ). 



4 Beyond the Glauber Approach. 

In this section we discuss the corrections to the Glauber approach (the Mueller formula 
of Eq.(66) ) as well as the way to construct a complete theory for deep inelastic scattering 
off a nucleus. 



4.1 The second iteration of the Mueller formula. 

To understand how big could be the corrections to the Glauber approach we calculate 
the second iteration of the Mueller formula of Eq.(66). As has been discussed, Eq.(66) 
describes the rescattering of the fastest gluon ( gluon - gluon pair ) during the passage 
through a nucleus ( see Figs. 1 and 2 ). In the second iteration we take into account also 
the rescattering of the next to the fastest gluon. This is a well defined task due to the 
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strong ordering in the parton fractions of energy in the parton cascade in leading ln{l/x) 
approximation of pQCD that we are deahng with. Namely: 



Xb < Xn < ... < Xi < 1 ] (89) 

where 1 corresponds to the fastest parton in the cascade. 

Therefore, in the second interaction we include the rescatterings of the gluons with 
the energy fraction 1 and Xi ( see Fig. 5 ). Doing the first iteration we insert in Eq.(66) 
G]\r{x,Q^) = G%^^{x,Q^). For the second iteration we calculate the gluon structure 
function using Eq.(66) substituting 

xG,, = "^^^^^'^'^ - xG^^^(x,Q2) ; (90) 

where xG\ is the result of the first iteration of Eq.(66) that has been discussed in details 
in section 3. 



Fig. 17 shows the need to subtract xG%^^ in eq. (|90|) making the second iteration. 
Indeed, in the second iteration we take into account the rescattering of gluon 1' - gluon 
2' pair off a nucleus. We picture in Fig. 16 the first term of such iteration in which G1/G2' 
pair has no rescatterings. It is obvious that it has been taken into account in our first 
iteration, so we have to subtract it to avoid a double counting. 



One can see in Figs.|T8| that the second iteration gives a big effect and changes crucially 
Ri, < 7 >, and < u >. The most remarkable feature is the crucial change of the value 
of the effective power a;(Q^) for the "Pomeron" intercept which tends to zero at HERA 
kinematic region, making possible the matching with "soft" high energy phenomenology. 
It is also very instructive to see how the second iteration makes more pronounced all 
properties of the behavior of the anomalous dimension ( < 7 >) that we have discussed. 
The main conclusions which we can make from Figs.|T8] are: (i) the second iteration gives 
a sizable contribution in the region x < 10^^ and for x < 10~^ it becomes of the 
order of the first iteration; (ii) for x < 10^^ we have to calculate the next iteration. It 
means that for such small x we have to develop a different technique to take into account 
rescatterings of all the partons in the parton cascade which will be more efficient than the 
simple iteration procedure for Eq.(66). However, let us first understand why the second 
iteration becomes essential to establish small parameters that enter to our problem. 



4.2 Parameters of the pQCD approach. 

As has been discussed, we use the GLAP evolution equations for gluon structure function 
in the region of small x. It means, that we sum the Feynman diagrams in pQCD using 
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Figure 18: Second iteration calculations for Ri, < uj >, and < 7 > for Ca and Au. 
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the following set of parameters: 



1 I 

as -C 1 ; as In- < 1 ; as In —2 < 1 ; as In — ^ In - ^ 1 . (91) 



The idea of the theoretical approach of rescattering that has been formulated in the GLR 
paper is to introduce a new parameter |^: 

- ^^"^"^a:G(x,g^) (92) 



2 



2Q2 

and sum all Feynman diagrams using the set of eq. (|9ll) and k as parameters of the 
problem, neglecting all contributions of the order of: as, as asln(l/x), asln(l/x) k, 
asln((5^/Qo) c>^s ln(Q^/Qo) ^- should be stressed that Mueller formula gives a so- 
lution for such approach. Indeed, Eq.(66) depends only on k absorbing all 
(as \n.{Q'^ IQl) ln(l/x) )" contributions in xG^x^Q"^). However, it is not a complete solu- 
tion. To illustrate this point let us compare the value of the second term of the expansion 
of Eq.(62) with respect to cr(r^) with the first correction due to the second iteration in 
the first term of such an expansion. In other words we wish to compare the values of the 
diagrams in Fig. 18a and Fig. 18b. The contribution of the diagram of Fig. 18a is equal: 

AxG{x,Q'){Ftg.l8a) = f^, J ^ J dQ" K^ix',^) ^ (93) 

where x' and Q'^ are the fraction of energy and the virtuality of gluon 1 in Fig. 18a. 
The diagram of Fig. 18b contains one more gluon and its contribution is: 

AxG{x,Q^)iFtg.l8b) = -0- / ^ ^ / — / ' ' 



71^2 IT J x' Q'^ J x" J ^ ^ ' 4 ' 
oc ln(l/x) IniQ^Ql) AxG{x,Q^)iFig.l8a) , (94) 

71 

where x' {x") and Q'^ ( Q"^) are the fraction of energy and the virtuality of gluon 1 
(1') respectively in Fig. 18b. Therefore, eq. (p^ ) gives the contribution which is of the 
order of eq. (|93D in the kinematic region where the set of parameters of eq. (|9lD holds. It 
means also that we need to sum all diagrams of Fig. 18b type to obtain the full answer. 
In the diagram of Fig. 18b not only one but many gluons can be emitted. Such emission 
leads to so called "triple ladder" interaction, pictured in Fig. 18c ( see ref. |]ll ). This 
diagram is the first from so called "fan" diagrams of Fig. 18d. To sum them all we 

* In the GLR paper the notation for k was W, but in this paper we use k to avoid a misunderstanding 
since, in DIS, W is the energy of interaction. 
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Figure 19: Corrections to the Glauber approach. 
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can neglect the third term in Eq.(65) and treat the remained terms as an equation for 
xGlXjQ"^). It is easy to recognize that we obtain the GLR equation Generally 
speaking the GLR equation sums the most important diagrams in the kinematic region 
where as ln(l/x) Ib.IQ'^/Qq) 3> 1 and k < 1. Using the Mueller formula we can give 
more precise estimates for the kinematic region where we can trust the GLR equation. 
Indeed, in Fig. 19 we plot the ratio 



R 



c 



7r22 



xGAix,Q^){Eq.{75)) - AxG%^''{x,Q^) ' 



(95) 



If Rc = 1, all the SC can be evaluated within good accuracy by the second term in 
Eq.(65). 
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Figure 20: The ratio Rq as a function of ln{l/x) and Q^. 



From Fig.^ one can see that for Ca we can safely restrict by the second term in the 
master equation and use the GLR equation to take into account the interaction of all 
parton in the parton cascade even for low values of in the HERA kinematic region 
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{x > 10 ^). However, for the nuclei we have to develop more general procedure for the 
interaction of the Mueller formula than the GLR equation for x < 10~^. 

We need to make some very important remarks, concerning the whole approach based 
on the Glauber - type shadowing corrections. It has been proven [^] that keeping all 
parameters of eq. ( ]90[) and eq. (pID and summing all Glauber - type interactions, as the 
Mueller formula does, is not enough. It turns out that the interaction between partons 
from different parton cascades that interact with the different nucleons are important. 
Fig. 20a shows the first interaction of such a type that has to be taken into account. This 
diagram should be compared with Fig. 20b which shows the interaction included in the 
Mueller formula as well as in the GLR equation. Fortunately, these new contributions are 
proportional to and we will neglect them, considering Nc is large enough. The general 

procedure how to sum all corrections of the order of at least for the GLR equation has 
been developed in ref. |^0[. In spite of the small parameter the parton selfinteraction 



can be essential for the case of a nuclear target but we postpone the detailed discussion 
of this problem to future publication. 




Figure 21: The Glauber approach ( Fig. 20b) and corrections to the Glauber approach due 
to interaction between partons from different parton cascades. 
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4.3 The GLR equation for nucleus. 



In this subsection we are going to discuss the GLR equation as a way of taking into 
account the interaction of all partons in the parton cascade with the target in spite of the 
criticism of the previous section. Indeed, let us go back to the discussion of the behavior 
of the average anomalous dimension in the Mueller formula and in the first iteration of 
the MF ( see Figs. 14 and 19). The general feature of the both iterations is the fact 
that the resulting value of the anomalous dimension ( < 7 > ) turns out to < 7 > < ^ 
for > IGeV"^. In this case the high order terms in the expansion of the MF are 
concentrated at large distances while the second one can be divergent at small distances 
if < 7 > — > |. Therefore we can rewrite the Mueller formula in the form: 



4 /■! fix' 

xGAix^Q') = ^ —I dhi (96) 



TT^ Jx X' Jo 



iff- 



irrf irrf 

Qq 



where Qo we choose of the order of 1 GeV"^. Recalling that < 7 > < 1 we can expand the 
first integral in eq. (^) and neglect all high contribution except the second term which 
can be important in the region where < 7 > |. Differentiating with respect to InQ^ 
and ln(l/a;) one obtain the GLR equation, namely 

S^.G.(..g^) _ iV.«s,^^(,_Q.) _ ^ {,:GA.,Q-)Y . (97) 



ain(l/x)91nQ2 ' ^ 

The second integral in eq. (^) gives the initial condition for the GRL equation. We 
want to draw your attention to two important outcomes from this simple consideration. 
First, the initial condition should be set only at sufficiently large value of Q^, e.g., at 
Q"^ = Qo — ^GeV"^. Second, we cannot use the MF to calculate this initial condition for 
the nucleus using the nucleon gluon structure function, since the corrections to MF for 
< IGeV"^ is large. Therefore, the solution of the GLR equation is only reliable in a 
kinematic region where it does not depend on any initial distribution. The way out of 
this shortcoming is only to use the direct information on the gluon structure function in a 
nucleus. Therefore, in such approach the main advantage of the Glauber formula is lost: 
the possibility to calculate the nucleus structure function from the nucleon one. This is 
the reason to develop a more general approach in the next subsection. 
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4.4 The generalization of the Glauber Approach. 



We suggest the following way to take into account the interaction of all partons in a parton 
cascade with the target. Let us differentiate the Mueller formula over y = ln(l/x) and 
i^HQVQD- It gives: 

. 14^3! { C + in. + } . (98) 

oyo^ 

Rewriting eq. (^) in terms of k given by 

Ncas 



xGa{x,Q^) (99) 



we obtain: 



o + — 7^ = {C + InK + Ei{K{y,0) I = FW ■ 100 

oyot, oy 71 

Now, let us consider the expression of eq. ( |100| ) as the equation for n. This equation 
sums all contributions of the order (as 1/0" absorbing them in xGA{y,0> ^^^^ all 



contributions of the order of k". For N^. ^ oo eq. ( p.00| ) gives the complete solution to 
our problem. The great advantage of this equation in comparison with the GLR one is 
the fact that it describes the region of large k and provides the correct matching both 
with the GLR equation and with the Glauber ( Mueller ) formula in the kinematic region 
where asy^ < 1. 

Eq. ( |1(JU| ) is the second order differential equation in partial derivatives and we need 
two initial ( boundary ) conditions to specify the solution. The first one is obvious, 
namely, at fixed y and Q"^ ^ oo 



The second one we can fix in the following way: at x = xq {y = yo) which is small, namely, 
in the kinematic region where asyC, < 1 

K ^ Kin = i^?t,2 xGa{x,Q'^) , (101) 

where xGa is given by the Mueller formula ( see Eq.(66) ). Practically, we can take 
Xo = 10~^, because corrections to the MF are small at this value of x = xq- 
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4.5 The solution to the generahzed evolution equation. 



4.5.1 The asymptotic solution 



First observation is tlie fact tliat eq. (|100|) lias a solution which depends only on y. Indeed, 

ition of the following equation: 

F{Kasymp) . (102) 



one can check that n = Kasympiy) is the solution of the following equation: 



^asymp 



dy 

The solution to the above equation is: 

sympiy) (if<i 



y-yo- (103) 



It is easy to find the behavior of the solution to eq. ( |103| ) at large value of y since 
F{k) asluK at large k ( as = "^as )• It gives 

i^asymp "sZ/ ln(as2/) "sZ/ > 1 ■ (104) 

At small value of y, F{n) — > aa,n and we have: 



^asymp *■ ^asymp 



{y = yo) e^^^y-y'^ . (105) 



The solution is given in Fig.^for as = 1/4 in the whole region of y for different nuclei in 
comparison with our calculations based on the MP. We chose the value of Hasympiy = yo) 
from eq. ( |101| ). We claim this solution is the asymptotic solution to eq. ( [IUCI| ) and will 
argue on this point a bit later. The calculations in Glauber approach for nucleon overshoot 
the asymptotic solution at large values of in HERA kinematic region ( at a; < 10~^ 
). However, at small values of the Glauber approach leads to stronger SC than the 
asymptotic solution. For nuclei the SC incorporated in the asymptotic solution turn out 
to be much stronger than the SC in the Glauber approach for any > 1 GeV^ at 
X > 10~^. In this kinematic region the solution of eq. (|10CI|) is drastically different from 
the Glauber one. 



A general conclusion for Fig.|22| is very simple: the amount of shadowing which was 
taken into account in the MF is not enough , at least for the gluon structure function in 
nuclei at x < 10^^ and we have to solve eq. ( p.OCI| ) to obtain the correct behavior of the 
gluon structure function for nuclei. 

Now, we would like to show that solution eq. QlOOD is the asymptotic solution of the 
new evolution equation. In order to check this let us try to find the solution to eq. ( |10(J| ) 



49 



300.0 
xGJA 



200.0 



100.0 



0.0 



Glauber Sol. Q' 


= 1.0 Gev' 




=10. Gev' 


Asymp. Sol 


= 1 . Gov' 




=10. Gev' 



Nucleon 




15.0 20.0 
y=ln(1/x) 



300.0 
xGJA 



200.0 



100.0 



Ca 




15.0 20.0 
y=ln(1/x) 




15.0 20.0 
y=ln(1/x) 



Figure 22: The Glauber approach and asymptotic solution for different nuclei. 
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in the form k = Kasymp + Ak, anticipating that An is small. If it is so, the following 
linear equation can be written for An: 

d'An{y,0 dAKjy,^ dF{K) 

dyd^ + = . (106) 

The general solution to eq. ( |106| ) one can find going to the Mellin transform in respect to 
^, namely: 

AK{y,0 = [ ^A«:(t/,z/)e^(«-«o), (107) 
where contour C is taken to the right of all singularities in u for function AK{y, u). 
For function AK{y, v) we have the equation: 

+ 1) =F\y)AK{y,v) , (108) 

where we denote F'{y) = ^^^^ \ K=Kasympiy)- Solving eq. ( |108| ) we obtain the answer: 

^<y,0 = I T^A'^M exp[z/(e -eo) + I' F\y')dy'] = (109) 

J c 2,17: y + i Jyo 

= I 7r~AK(z/) exp[z/(^ - ^o) + \VL{F{Kasymp{y)) , 

J C 2,11: Z/ + i 

where function Ak.{v) should be find from the initial condition at = ^o? namely Ak(?/, = 
0) = 0. 

To satisfy the initial condition we will find function A/t(z/) from the equation: 

^^asymp{y.i = io) = / -— A/s:(z/) exp[ — — / F'{y')dy']. (110) 
JC 2ni I' + i Jyo 

In doing so we obtain Ak(j/,^): 

r dv 1 ry 

^<y.i) = / 77-Ak(z/) exp[^ / F'{y')dy']{eMH^-^o)) - 1}, (HI) 

Jc 2t7T U + 1 Jyo 

which satisfies all our requirements. 
We claim that function: 

r dv 1 rv 

<y,0 = / IT-A'^M exp[^ / F'{y')dy'] exp( - ^o) ) (112) 
Jc 2nx z/ + i Jyo 
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is the approximate solution to the nonhnear equation ( |10CI| ). 



Let us check this point considering very large values of ^ and y. In the region of large 
y, IyoF'{y')dy' as ln(as|/) and Ak(z/) oc ^ which reproduces Ka^y^p ^ asl/ at large 
y and ^ = ^o- Making use of the saddle point approximation one can see that 

AKiy,0 ^ eV€"sin(f>a™si/) « (113) 

To give a solution able to describe the experimental data we have to adjust the behavior 
of i^asymp at Small y {y < yo) with the available parameterization of the gluon structure 
function, in particular with the GRV parameterization which we have been using through 
this paper as a standard one. We have not done this job in this paper and we are going 
to publish the result of this calculation elsewhere. However, to estimate the effect of the 
SC which follow from eq. (|100|) we solve this equation using the semiclassical approach. 



4.5.2 Semiclassical Approach 

Here we solve eq. ( P-OOp using the semiclassical approach, adjusted to the solution of the 
nonlinear equation of eq. ( |10(]| )-type in refs. ^ . For simplicity, we assume that as 

is fixed. 

In the semiclasical approach we are looking for the solution of eq. ( |100| ) in the form 

K = (114) 

where 5* is a function with partial derivatives: ^ = lu and fl = 7 which are smooth 
function of y and ^. It means that 

d^S dS dS 

<T^-^ = ^7 (115) 



d^dy dy 

Using eq.( |115| ), one can easily rewrite eq. ( |100| ) in the form 

dS dS dS 



or 



c^(7 + l) = $(5) (117) 
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We are going to use the method of characteristics ( see, for example, ref. |^). For 
equation in the form 



F{^,y,S,^,u;) = 



'1181 



we can introduce the set of characteristic hues {^(t),y(t), S(t),u!(t),'j(t)), functions of the 
variable t, which satisfy the following equations 



„^ dy 

— = Fy] — = F^; — = ^Fy + uF^ ; 



dt 

dj 



-{F^ + jFs 



dS_ 

Itt 
duo 

~dt 



-{Fy + LuFs 



(119) 



where F^ = ■^F{^,y, S,'j,uj) etc. Eg . (11191) looks as follows for the case of eq (|117|) 



df dy dS , , 

^/ ^1 



(120) 



where $5 = H. For practical purpose it is better to rewrite the set of equations ( |120|) in 
the form 



d^ 



dS 



uo- 



(27 + 1) rf7 



dy 7 + 1 dy 7 + 1 ' c/y 

Using eq.( |117|) , eq.( |121|) can be rewritten 



7 



'7 + 1 



:i2r 



di $(5) dS 27 + 1 



(i7 



dy (7 + 1)2' + ^ dy S + 1 

The initial condition for this set of equations we derive from eq.( |101j ), namely 



7 



(122) 



7o 



So = IriKiniyo^io) 
dlnKin{yo,0 



(123) 



Let us discuss the main properties of the solution before numerical calculations. The 
first observation is that $5 < and $(5) > for all values of S. From the second 



equation of the set ( |122| ) we see that S decreases along all trajectories with 7 < —1/2 
and increases for trajectories with 7 > —1/2. Thus, it is useful to study the qualitative 
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behavior of the trajectories for two different regions of the initial condition. Namely, for 
7o < -1/2 and for 70 > -1/2. 

From the third equation of the set (|122|) , we notice that ^ > for all 7 < 0. It 
means that 7 grows with y starting from 79. However for 7 > 0, ^ < and 7 starts to 
fall down. In both cases, for 70 > —1/2, 5* goes to infinity as y grows, and $5 and the 
derivative ^ go to zero. Thus, we can conclude that 7 ^ as y ^ 00. 

It is useful to study closer what is happening at small 7. 

From equations ( |122| ), we can write 

dS_ _ 27 + 1 ^S) 
d^"7(7 + l)<f'5(5) 

dln<^ 27 + 1 



(124) 



or 



(i7 7(7 + 1) 



(125) 



which has the solution 



For small I7I the solution of eq. (|126|) is correlated with large values of 5", since ^{S) 
has maximum at S* = 5*0 and decreases at S" > 5*0. On the other hand, at large S", 
$5 — >• — $(S'). Substituting this relation in the third equation of (|122|) , we obtain 

d(7/7o) $(5o) 7' 



dy 70 + 1 % 



2 ' 



(127) 



which solution is 



1 +l = _i(M(,-,„) (128) 



or 



7/70 7o + 1 



7 



70 ^(y-yo) + i 



(129) 



Therefore, we see that 7 approaches 7 = at large j/ for 70 > —1/2, either for 7 positive 
or negative. 
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In order to find the solution of wliole set of equations ( |122| ) at large values of y we use 
the fact that 7^1. Indeed, neglecting 7 in comparison with 1, eq's. (|122|) are reduced to 



•S* — — ^ — ^0 

§^ = ^{S) . (130) 
dy 

Rewriting the second equation in terms of the function k = e"^ we have 

l^fW. (131) 

The solution of this equation is k = t^asymp- The first equation gives the equation for 
trajectories and ai y ^ yQ we have 

^ - ^0 = /n (^ ^"" ^ Inasiy - yo) + Mnasiy - yo) (132) 



These trajectories are the same as the trajectories of eq. (|106|) . The simplest way to see 



this is just to find the saddle point in the solution of eq.(|109D. 



Therefore the qualitative picture of the trajectories looks as follows. At small values 
oi y — yo we can start from initial condition in which e"^" ^ 1. In this case $5^^ oc — e'^", 
g^o ^ I g^j^^ ry remains close to 70 in the large interval of y — yo. For these values of y — yo, 
S grows as a function of y and this grows leads $5 to approach zero. This behavior reflects 
in the decrease of 7 versus y — yo. Finally, at very large values of y — yo we approach the 
asymptotic solution for large region of ^. 

Now we will study the qualitative behavior of the trajectories of eq. (|122|) for 7 < —1/2. 
As we already know, S decreases along all trajectories with 7 < —1/2. As S goes to 
negative values, $(5*) goes to and $5 goes to zero from negative values. It means 

that d'j/dy > but tends to zero as y grows. Thus, 7 increases and tends to a constant. 
As S decreases, the value of k goes to zero. This behavior is enhanced for 70 closer to —1. 
We will show below that this solution approaches the solution of the GLAP equation in 
DLA for 7o < -1/2 

Let us recall that in the GLAP equation $5 = and $(S') = Thus, the set of 

eq's(|122[) can be rewritten in the form: 



d^ Ncas _ 

(133) 



dy 7r(7o + 1)2 
dS Ncas 270 + 1 



dy vr (70 + 1)2 
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where 70 is the initial value of the 7 which does not change along the trajectory. The 
solution of eq (|133D is 



7r(7o + ly 

S-So= .^^"^^., (270 + l)iy- yo) ; (134) 



or 



7o + 1 



I Ncas y- yo 



S-So = -{^- eo) + 2 J^(^ - ^o){y - yo) ■ (135) 



The above expression is the solution of DLA equation for k. From eq.( p.34| ) we see that k 
in DLA approximation goes to zero for 70 < —1/2 

In the region of the double log approximation we consider 

Ncas 



(e-eo)(y-yo)~l (136) 



while ^^(^ - ^0) < 1 and ^^{y - yo) < 1. Therefore 

70 + 1 ^ < 1 (137) 

We set the initial condition y = yo = 4.6 {xb = 10~^), where the shadowing correction is 
not big and the evolution starts from 7 < 0. In this case d'j/dy > and the value of 7 
increases . At the same time dS/dy < and S decreases if 70 < —1/2. With the decrease 
of S, the value of $5 becomes smaller and after short evolution the trajectories of the 
nonlinear equation start to approach the trajectories of the GLAP equations. We face 
this situation for any trajectory with 70 close to -1. If the value of 70 is smaller than — | 
but the value of 5*0 is sufficiently big, the decrease of S due to evolution cannot provide a 
small value for $'(5") and 7 increases until its value becomes bigger than — | at some value 
of y = yc- In this case for y > y^ the trajectories behave as in the case with 70 > — 1. 
For 7o > —1/2, the picture changes crucially. In this case, dS/dy > , d'j/dy > and 
both increase. Such trajectories go apart from the trajectories of the GLAP equation and 
nonlinear effects play more and more important role with increasing y. These trajectories 
approach the asymptotic solution very quickly. 
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For the numerical solution we use the 4th order Runge - Kutta method to solve our set 



of equations with the initial distributions of eq. (|123|) . The result of the solution is given 
in Figs.^ and In these figures we plot the bunch of the trajectories with different 
initial conditions. For the nucleon ( Fig.^ ) we show also the dependence of 7 along 



these trajectories. One can notice that the trajectories behave in the way which we have 
discussed in our qualitative analysis. It is interesting to notice that the trajectories, which 
are different from the trajectories of the GLAP evolution equations, start ai y = i/q = 4.6 
with the values of between 0.5GeV^ and 2.5GeV^ for a nucleon. It means that, guessing 
which is the boundary condition at = Qq = 2.5GeV^, we can hope that the linear 
evolution equations ( the GLAP equations) will describe the evolution of the deep inelastic 
structure function in the limited but sufficiently wide range of Q^. In other words, we can 
repeat the trick done in subsection 4.3 deriving the GLR evolution equation for nucleus. 



In Figs. p3| and El we plot also the lines with definite value of the ratio R 



xG(x,Q^)(qeneralized equation) /i • j i t \ rpi t • j.i j. ^- j. ^ 

— xG{x Q'^){GLAP) (horizontal hues), ihese hues give the way to estimate how 

big are the oC. One can see that they are rather big. 
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Figure 23: The trajectories and contour plot for the solution of the generalized evolution 

Pmmiin'n fnr N R — a:G(x,Q^)(generafi2:ed equation) 

equauon jot i\. n - xG{x,q^){glap) 



We have discussed only the solution with fixed couphng constant which we put equal 
to as = 0.25 in the numerical calculation. The problem how to solve the equation with 
running coupling constant is still open. The GLR equation is the limited case of the 
general one when we consider only two first terms in the expansion of F{k,) with respect 
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Figure 24: The trajectories and contour plot for the solution of the generalized evolution 

i- e J A D xG{x,Q^)(qeneralized equation) 

equation for La and Au. K = — xG(x q^){glap) 

to K. In Fig.^ we picture the trajectories and contour plot for the GLR equation. One 
can see that the GLR equation gives stronger SC that the generahzed evolution equation. 
However, for nucleon the difference becomes sizable only at very small values of x out of 
the HERA kinematic region. 

We would like also to draw your attention to the fact that our asymptotic solution 
turns out to be quite different from the GLR one. The GLR solution in the region of 
very small x leads to saturation of the gluon density p8| , ^ Saturation means that 



K, tends to a constant in the region of small x. The solutions of eq. (|100|) approach the 
asymptotic solution at a; — > 0, which does not depend on Q^, but exhibits sufficiently 
strong dependence of k on x ( see Fig.|22| ). 



5 Conclusions. 

1. The Glauber approach to the gluon structure function in a nucleus, suggested by 
Mueller in Ref. , has been developed and studied in detail. Using the GRV parameter- 



ization for the gluon structure function in a nucleon, we calculated the value as well as 
energy and dependence of the gluon structure function in a nucleus.lt is shown that 
the shadowing corrections are important in the region of small x and changed crucially 
the value and anomalous dimension of the nuclear structure function, unlike the nucleon 
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one. The interesting observation is the fact that the average anomalous dimension for 
bigger that IGeV^ does not exceed | for nuclei. This fact makes the application of the 
GLAP evolution equations much more reliable in the case of nuclei than for a nucleon. 

2. The corrections to the Glauber approach have been investigated in a systematic 
way and it turned out that they are essential and should be taken into account. The 
QCD technique has been developed to go beyond the Glauber approach and the new 
evolution equation for nuclear structure function has been derived in perturbative QCD 
( see eq. (IHO) )• 



3. The new evolution equation was solved in the semiclassical approximation and the 
main qualitative properties of the solution were discussed. In particular, it was shown 
that this solution does not lead to the saturation of the gluon density but manifests a 
strong dependence of gluon density on x. We consider that this solution can provide a self- 
consistent interface between "soft" high energy phenomenology and "hard" QCD physics. 
Surprisingly, the solution of the GLR equation gives a very accurate approximation to 
the solution of the new evolution equation in the HERA kinematic region. 

Several problems remain beyond the scope of this paper. First, we described the parton 
cascade in the GLAP evolution at low x or in other words in double log approximation 
of perturbative QCD (DLA). At first sight it is a controversial assumption since the DLA 
works rather badly in the accessible region of x and should be replaced by the BFKL 
dynamics at very small values of x. We found the possible way out of this difficulty since 
the SC change the value of the anomalous dimension in such a way that it does not reach 
the value of | where the BFKL equation should be taken into account. It means that the 
SC come first and the BFKL evolution will never develop. We consider this result as a 
plausible explanation why the BFKL evolution has not been seen yet in nucleon data at 
HERA. The experiments with nuclei will certainly shed light on this problem. 

The second problem that we have not discussed in this paper is the corrections to the 
Glauber approach due to the selfinteraction of the partons belonging to different branches 
of the parton cascade. We suppose to do this in further publications using the general 



approach proposed in ref. 40 



A certain limitation of all calculations in the paper is the fact that we used only the 
GRV parameterization of the gluon structure function in the nucleon. We would like only 
to recall that the goal of the paper is not to provide a reliable prediction for an experiment 
but rather to study the mechanism and amount of the shadowing corrections. Therefore, 
the GRV parameterization was a tool for our theoretical experiment which is based on 
the available experimental data and on the GLAP evolutions. Nevertheless, we are going 
to make full calculations using all parameterizations on the market in the nearest future. 
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We also used the simplest assumption for the nucleon density in a nucleus, namely,the 
Gaussian one. The calculation with more general parameterization of the nucleon density 
will come out soon. 

We hope that our paper will convince the reader that the SC are essential for the 
gluon density in a nucleus, and that the Glauber approach, in spite of the fact that it is 
widely used, is not enough both from the phenomenological and theoretical point of view. 
Fortunately, QCD gives us the first example how to treat the corrections to the Glauber 
approach theoretically. 
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